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We study the orbital angular momentum structure of the quarks inside the proton. By employing 
■ the hght-cone diquark model and the overlap representation formalism, we calculate the chiral-even 
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generalized parton distribution functions (GPDs) Hq(x,^, A ), Hg{x,^, A ) and Eg(x,^,A ) at zero 
^ , skewedness for q — u and d quarks. In our model and have opposite sign with similar size. 



Those GPDs are applied to calculate the orbital angular momentum (OAM) distributions, showing 
that Lu{x) is positive, while Ljix) is consistent with zero compared with Lu{x). We introduce the 
impact parameter dependence of the quark OAM distribution. It describes the position space dis- 
tribution of the quark orbital angular momentum at given x. We found that the impact parameter 
dependence of the quark OAM distribution is axiaUy symmetric in the light-cone diquark model. 
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■ I. INTRODUCTION 

oo 

CN \ understanding the spin structure of the nucleon is one of the most important challenges in hadron physics. 

Q ■ The naive picture that the nucleon spin is provided totally by the spin of its three valence quark was proved 

o ■ 



X 



to be wrong by the experimental measurements. The EMC result 11|] indicates that a large fraction of 
the nucleon spin is carried by other sources of angular momentum. There have been many attempts to 
explain the EMC result from the fundamental theory. Besides the angular momentum of the gluon, the 
quark orbital angular momentum (OAM) |2] is believed to provide a substantial part of the nucleon spin. 
In the last two decades the theoretical description of the quark OAM distribution has been established jsl- 
Q]. It has been shown by Ji that the quark angular momentum can be separated into |4] the usual quark 
helicity and a gauge-invariant orbital contributions Lq. One of the advantage of this decomposition is that 
Lq is related to generalized parton distributions (GPDs) Is MlSl] . the experimental observables that enter 



1411 and meson 



the descriptions of hard exclusive processes, such as deeply virtual Compton processes 
exclusive production 

Moreover, recently it has been found that the quark OAM plays an essential role through spin-orbit 
correlations in some novel phenomena that appear in the physics of single spin asymmetries, among which 



a particular transverse momentum distribution (TMD) \V1 , 



180 — Sivers function 119L |20|] — has attracted 



2 



a lot of interest, since it is an essential piece in our understanding of the single spin asymmetries (SSA) 
observed in semi-inclusive deeply inelastic scattering (SIDIS). These SSAs have been measured recently 



by both the HERMES Ul 



m and COMPASS ||23, 



there is a quantitative relation 125 



2411 Collaborations. An interesting observation is that 



2711 between the Sivers function f^^ and the GPD E'^, although it is 
obtained in a model dependent way, suggesting that similar underlying physics plays a role for nonzero 
ft~r!! and E"^. Similar relations have been obtained between Boer-Muldes functions and chiral-odd quark 
GPDs II28L 12911 ■ A complete study on the relations between the GPDs and TMDs has been presented in 1 30 ], 



which becomes more transparent through the conception of general parton correlation functions 



The relations between GPDs and TMDs are more intuitive 133 



111, 



position (impact parameter) space |135 



3411 if we interpret GPDs in the transverse 



3811 ■ Of particular interest is the case of zero skewedness = 0), 



where a density interpretation of GPDs in impact parameter space may be obtained 13511 . In particular this 
interpretation allows one to study a three-dimensional picture of the nucleon. 

In this paper, we study the orbital angular momentum structure of the quarks inside the proton in a 
light-cone diquark model. In this model the light-cone wave function of the proton can be obtained. It 



3,Q. 



is then convenient to express the physical observables in the overlap representation formalism ||3S 
We calculate the chiral-even generalized parton distribution functions (GPDs) Hq(x, ^, A^), Hq(x, ^, A^) and 
Eq{x, ^, A^) at the zero skewedness forq = u and d. We found that and E^i have opposite sign with 
similar size in this model. The GPDs are applied to calculate the quark 0AM distributions, showing that 
L„(x) is positive, while Ld{x) is consistent with zero compared with L„(x), and the net 0AM of the u and 
d quarks is positive. We also introduce the impact parameter dependence of quark 0AM distribution. It 
describes the position space distribution of the quark 0AM at given x. We found that the impact parameter 
dependence of quark 0AM distribution is axially symmetric in the light-cone diquark model. 

The manuscript is organized in the following way. In Section. II we review the GPDs and their con- 
nections with quark orbital angular momentum, in Section III, we present the calculation of chiral-even 
GPDs from the light-cone diquark model, by applying the overlap representation formalism. We also show 
the calculation of the quark 0AM in the same approach. In Section IV we introduce the impact parameter 
dependence of quark 0AM distribution and present results of the position space distribution for orbiting u 
quark, from the light-cone diquark model. We summarize our paper in Section V. 
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n. SYSTEMATICS OF GENERALIZED PARTON DISTRIBUTIONS AND THE ORBITAL ANGULAR 

MOMENTUM 

GPDs are introduced to describe the exclusive process in which the momenta of the incoming and 
outgoing nucleon in the symmetric frame are given by 

p = p+'J^, p' = p-^A, (1) 

and satisfy = p'^ = M^, with M denoting the nucleon mass. The GPDs depend on the following variables 

x=—, ^ = , t = A^, (2) 

p+ ^ 2P+ 

where the light-cone coordinates are defined by 

a- = {J'±a^), dT = {a\a^) (3) 

for a generic 4-vector a. In a physical process the so-called skewness ^ and the momentum transfer t to the 
nucleon are fixed by the external kinematics, whereas x is typically an integration variable. 

The chiral-even GPDs Hq, Eg and Hq, Eq for quarks are defined through matrix elements of the bilinear 
vector and axial vector currents on the light-cone: 



y+=0,>'x=0 



J S/r 

— U{p') I y^H{x, ^, t) + -J^E{x, ^, t) J U{p) , 

rdy_ 
J 

^ U{p') I y^ysHix, ^ E{x, f , t) J U{p) 



IP-" 

^^^/.vPV/2^^'|^(0)rV5<A(3')l/'> _ „ (5) 



2P-^ 

An important implication of GPDs is that they are related to the 0AM (0AM ) of the quark, which is 
expected to provide essential contribution to the total spin of the nucleon. Here we follow the decomposition 
of the nucleon spin introduced by Ji ||4|] : 

q q 

where Aq, Uq and 7| denote the quark spin, quark 0AM and gluon angular momentum, which comes from 
the expectation value of the operator 

= ^ 2 '^I^''^'? + Z 'I'l^^ '^^'"^^ + [fx{ExB)f, (7) 

1 9 
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FIG. 1: The generalized parton distributions Hu{x,Q, Aj) and Hd(x,Q, Aj) for the proton in the Hght-cone diquark 
model as functions of x for different values of At- 



Note that in literature |l3L|g,|4l|] there are some other ways to decompose the nucleon spin. The advantage of 
the decomposition of to Aq and L| in ^ is that it ensures the gauge invariance of the operators. There has 
been also discussion that wether the gluon angular momentum can be further decomposed gauge-invariantly. 
In this work we will not consider the gluon contribution. 

The quark 0AM distribution Lq(x) can then be defined as the expectation value of operator 



between the proton state \PS}: 



Ll{x) = {PS\dL\PS) 



(8) 



(9) 



The quark 0AM distribution can be obtained from jA , 



42|] 



L^(x) = Ux [Hgix, 0, 0) + Egix, 0, 0)] - H^ix, 0, 0)) , 



(10) 



where Hq{x, 0, 0), Hq{x, 0, 0) and Eq{x, 0, 0) are the forward limits of GPDs. Furthermore, the former two 
are the unpolarized and hehcity distributions for the nucleon, respectively. 



q{x) = Hq(x, 0, 0), Aq{x) - Hq(x, 0, 0), 
and Eq{x, 0, 0) is related to the anomalous magnetic momentum of the nucleon in the following way: 

I dxEq{x, 0, 0) = Kq, 

Jo 

where Kq is the contribution of quark flavor q to the nucleon anomalous magnetic momentum. 



(11) 



(12) 
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FIG. 2: The generalized parton distributions Hu{x, 0, A^) and Hd(x, 0, A^) for the proton in the light-cone diquark 
model as functions of x for different values of Ar. 



III. GPDS IN THE LIGHT-CONE DIQUARK MODEL FROM THE OVERLAP REPRESENTATION 

FORMALISM 

In this section we present the calculation of the GPDs in the light-cone diquark model from the over- 
lap representation formalism. The proton wave function with helicity tl, JJ. in the SU(6) quark-diquark 



model ll43l445l] in the instant form is written as 



^i^\qD) = J-^y\qVf'^ + -L^slqS)^'^, (13) 
where D = V,S denotes the vector diquark and scalar diquark, respectively. The 

- ^Vo{uu)d^'^ + 2V±i{uu)d^''^]; 
\qSf-^ = SiudW'K (14) 

The spin part of the light-cone wave function of the proton can be obtained from the instant form of the 
wave function by a Melosh rotation. For a spin- 5 particle, the Melosh transformations are known to be ||46|] 



X^j = co[{k^ + m,)x^p + k'^xl], (15) 

where ;k'7- andxF are instant and light-cone spinors respectively, to - ^Ik''' (k^ + m^^j ^ , k^-^ = k^ ± ik^, and 
niq is the quark mass. In this work, for simplicity we treat the diquark as a point particle. The scalar diquark 
does not transform, since it has zero spin. For the spin-1 vector diquark, the Melosh transformations are 
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given by 14711 



col [ V2 + Av) 4 + 2 ((^° + Ay) - k^k^) V] 
-^f2{k^y + Av)k^V-p'], 

+ V2 + Av) k'^V'j, + (k^ + Avf V^^ 



(16) 



Here, Ay denotes the mass of the diquark, Vt and Vp are the instant and light-cone spin-1 particle respec- 
tively, which are constructed within the Weinberg-Soper formalism |48]. 

After some algebra we arrive at the two body light-cone wavefunctions of the proton with 



1 



1 



1 



V2 V6 V3 

The scalar diquark component of the wavefunction for the proton has the form 



X ^l''^{x,kT,Sz)\xP'',kT,s,}, 



(17) 



(18) 



while the vector diquark component is expressed as 



d^kjdx 



Vx(l -x)167r3 
X i]jI\x, kr, k, s,)\xP*, kr, k, s,), 



=0,±l;i,=±^ 



(19) 




0.2 



J I I I I I I I L 



A.j,=0 GeV 
Aj,=0.4 GeV 
A.j,=0.8 GeV 
A.|,= 1.2GeV 
Aj,=1.6GeV 



0.4 0.6 

X 



0.8 



FIG. 3: The generalized parton distributions Ei,{x,0, A^) and Ed{x,0, A^) for the proton in the light-cone diquark 
model as functions of x for different values of Aj-. 
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which is the same for \u V)f and \d V)f. Here we denote s, and as the spin projections of the quark and 
the vector diquark. The forms of ij/^ {x, kj, Sy) and ifjy\x, kj, 1^, Sy) are given in the appendix. 

Now we calculate the chiral-even GPDs in the zero skewedness = 0) where t - -A^. In the overlap 



representation [391, |4(J|] H, E and // at ^ = can be expressed in a symmetric frame as ( in the domain 
< X < 1 and for « — > « transition): 



7^(x,0,-A|) 



n,Aj i=\ 



n,Ai i=\ 



















5(2) 




















( " 






r » ^ 










^(2) 





















H(x,0,-Al) = 2] / Y\^ign(A,)^^^l6n^l-f^xM^ 

n,Ai i=l V ;=i / w=i 



with 



(22) 



x[ = xi, k^^ = kri - (l - xi)— 
for the final struck quark, 

Xj — Xj, kjj — kji + ^'~2~ 

for the final (n-1) spectators. 



and 



- kji +(1 -x\)— 



yi = x\, /xi 

for the initial struck quark, 

Ji = Xi, hi 

for the initial (n-1) spectators. 



_ Ay 
f^Ti Xi 2 



From Eq. (|2TI) we see that non-zero Eq needs a spin flip between the initial and final proton wavefunc- 
tions. The same kind of overlap integration of light-front wavefunctions (with = ±1/2 in the initial and 
final states) also appears in the calculation 149(1 of Sivers functions, which indicates the presence of the 
quark 0AM . 

By employing the light-cone wavefunctions given in (fTTl) and the overlap representation formal- 



ism, we calculate the generahzed parton distribution functions (GPDs) Hq{x,0,-Ai), Hq{x,0,-\l) and 
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FIG. 4: The OAM distributions Lq{x) of u and d quarks inside the proton in the light-cone diquark model as functions 
of X. 



Eq{x,0, at zero skewedness for ^ = m and d quarks. The x-dependence of these GPDs at different 
values of Ay are given in Figs. [T]|2] and [3j respectively. 

From Fig. [3] one can see that £"„ and have opposite sign is positive and £^is negative) with 



similar size in our model. Since it has been shown that there is a quantitative relation 
the Sivers function /j"^^ and the GPD our ...coincdeswiU, recent ex.acUo„.Ua of the Sivers 



121121 

dm 



3011 between 



function from the Semi-inclusive deeply inelastic scattering data, which show the Sivers functions of u and 
d have opposite sign with similar size. 

Special attention should be paid to the limit of zero momentum transfer = 0, since in this limit 
the GPDs and Hq are simplified to the forward distribution q{x) and Aq{x). Also the quark OAM s 
are related in the way shown in (ITTI ). from which in principle one can calculates Lg{x) from the known 
chiral-even GPDs. 

By taking the GPDs in the forward limit, we calculate the OAM distributions of u and d quarks inside 
the proton, as shown in Fig. ID It can be seen that in our model L„(x) is positive, while L^ix) is consistent 
with zero compared with Lu{x), and the net OAM of the u and d quarks is positive. From Fig. [3] one can 
see that Ed is sizable. However, since Ed is negative, there is a cancelation between d{x), Ed{x) and A^(x). 
This leads to a small contribution of the d quark orbital angular momentum. We remind that there are 



Lattice QCD 153 



1155 



5411 . as well as phenomenological parametrizations and other model calculations of GPDs 



5911 . which are used to estimate the OAM of the quarks. 
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FIG. 5: The impact parameter distributions (scaled with a factor of i2n)^) xL^ix, bj) (left) and xL^ix, bj) (right) for 
the proton in the light-cone diquark model as functions of x for different values of b. 




FIG. 6: The profiles of the impact parameter distribution (scaled by a factor of (2n)^) xLu(x, bj) for the proton in the 
hght-cone diquark model as functions of Aj- for x - 0.3 (left) and x - 0.5 (right). 



IV. IMPACT PARAMETER DEPENDENCE OF ORBITAL ANGULAR MOMENTUM 



In this section we want to study the quark 0AM s in transverseposition (impact parameter) space. The 



GPDs in the impact parameter space have been studied in Refs. 1135 



371] ■ The most interesting case is the 



zero skewedness limit ^ = 0, in which a density interpretation of GPDs in the impact parameter space 



may be obtained i35h . Therefore studying GPDs in impact parameter space can provide a three-dimensional 
picture of the nucleon. In the following we restrict ourselves to the case ^ = 0. 

The impact parameter PDFs inside the nucleon can be obtained by sandwiching the parton correlator 



between nucleon states localized in transverse space 

q{x,bT) = {p\OT;s\df\x,bT)\p^,(iT;S}, 



(23) 
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where 



and the initial and final states in the transverse space defined as 135 , 

72 



31 



\P\br,S) = N r ^^>-'>-*-|p;5), 



{P*,bT;S \ =N' 



r (fp' . 

J {2nf 



(24) 

(25) 
(26) 



which characterize a nucleon with momentum at a transverse position bj and polarization specified by 
S. 

One of the interesting features of impact parameter dependent parton distributions is that they are Fourier 



transformations of GPDs i35h . For instance, The impact parameter dependence of unpolarized quark in the 



unpolarized nucleon can be obtained from 

q{x, bT)= \ — 

J (27r. 



) 



(27) 



here bj and are two conjugated parameters. 

Similarly the impact parameter dependence of quark helicity distribution in the longitudinal polarized 
nucleon is defined as 



Aq{xM) = {P\OT;s\df'"\x,bT)\P\OT;S), 



Which is the Fourier transformation of H^: 



q{x, bj) 



Cd^A 
"J (27:. 



|^.-'*-^-//,(x,0,-A|), 



(28) 



(29) 



We follow a similar approach by introducing the impact parameter dependence of quark 0AM br). 
It can be obtained from the expectation value of Ol, given in Eq. dSjl, between the position state IP""", Or): 



L^{x,bT) = {POt;S IOlIPOtIS) 
After a Fourier transformation on Lq{x, bj) one can anive at 

J d^br e'*^ ■ L\{x, br) - L^ix, - A^ ) 

The function L|(x, -A|) can be obtained by the GPDs at zero skewedness 

L^(x,-A|) = ^{x[h{x,0,-\1) + Eg{x,0,-Al)] 
- H{x,0,-Al)], 



(30) 



(31) 



(32) 
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and (fTTI ) is the forward limit of Lq{x, -A^). 

Therefore, if we know the GPDs Hq Hq and Eq, from (fTTI ) one can calculate the impact parameter 
dependence of the quark 0AM distribution by the Fourier transformation 

^e-'*-^^L^(x,A2). (33) 

The integration over impact parameter dependence of quark 0AM leads to 

J d'bTLq{x,bT) = Lq{x) (34) 

In Fig. 15] we shown the impact parameter distributions (scaled with a factor of {In)^) L^ix, br) (left) 
and Lrf(x, bj) (right) for the proton in the Ught-cone diquark model, as functions of x, for different values 
of b. In Fig. |6] we show the profiles of the impact parameter distributions L„(x, br) for the proton in the 
light-cone diquark model as functions of bj, for x = 0.3 and x = 0.5. It is shown that the impact parameter 
dependence of quark 0AM is axially symmetric. Also at large x the impact parameter distribution is peaked 
at small b. 



V. SUMMARY 



As a conclusion, we study the 0AM structure of the quarks inside the proton in a light-cone diquark 
model. In this model the light-cone wave function of the proton is known. It is then convenient to express 
the physical observables in the overlap representation formalism. We calculate the chiral-even generalized 
parton distribution functions (GPDs) Hq{x, f , A^), Hq{x, ^, A^) and Eq{x, ^, A^) at zero skewedness for 17 = m 
and d. We found that and E^ have opposite sign, with similar size in our model. The GPDs are applied 
to calculate the 0AM distributions, showing that Luix) is positive, while Ld{x) is consistent with zero 
compared with L„(x), and the net 0AM of the u and d quarks is positive. We also introduce the impact 
parameter dependence of quark 0AM distribution L(x, br) ■ It describes the position space distribution of 
the quark 0AM at given x. We found that the impact parameter dependence of quark 0AM distribution is 
axially symmetric in the light-cone diquark model. 
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Appendix A: light-cone wave functions in a diquark model 



The expressions for i/f^'^(x, kj, have the form 



il/l{x,kT,+-) - (ps{x,kT), 



(jj 



il/l{x,kT,-\) - -—4>s{x,kT), 

L to 



kr 



and 



ll/l{x,kT,+^) ^ —(f>s{x,kT), 



i//Ux,kT,—^) 



1, 
r 

2' 



Cl) 

(yfc+ + m) 
a> 



(psixykr). 



respectively. 

The expressions of i/fy^ix, kr, Iz, s^) can be expressed as 

-4>v{x, kr) 



ilfl(x,kT,+l,V = -V2^ 



OJOJ', 

+(kt + Av) 



{k^ + AvXk"- + m) 



k\ 



OJOjfr L 



-{k^ + Avfik^ + m) 



il,\{x,kT,%^) = 2 



(jJCO 



{kl + Av)k-^ - kl 



+ m) 



4(x,fcr,0,i) = 



-{k^ + Av)k^, 
(l>vix,kT) 



0J0J7 



2((4 + Av)k;, - 4) 



4(x,*r,-l,T) = V2 



2{klr + Av){k^ + m) 
0y(x, kr) 



k\ 



l/fl{x,kT,-hi) = 



(f>v(x, kr) 



(k^ + AvXk"- +m)-k\ 
l4\k^ + Av + k^ +m) 



k^. 



(Al) 



(A2) 



(A3) 



13 



and 



(At(x,A:r,+l,t) 



kj(ky + Av + k'^ +m) 



(kt, + Av){k^ + m) 



k\ 



ilj\{x,kT,OA) 



iijl{x,kT,Q,i) 



= 2 



(pvix, kr) 



({k^y + - kj) 



+{ky + Av){k'^ + m) 



k\ 



6v{x, kj) 



0)0), 



{k^y + Av)k^ - k\ 



+ m) 



{k^ + Av)kU, 



4(x,A:r,-l,T) = V2 



(pv{x, kr) 



OJO) 



V 



{kl + Av)k\ 



-{k^ + Avfik^ + m) 



iffl{x,kT,-l,i) = V2 



(pv{x, kr) 



OJO) 



V 



Hkv + ^vY 



(/c+ + Av){k'' + m) 
k', 



(A4) 



The momentum dependence of the wavefunctions in the above equations is described by (poix, kj) with 
the Gaussian form 



(^z)(Jc, kr) = Ad exp 



32 



where 



kl. + ml ki + aI 



9 _l_ "T ' 



(A5) 



(A6) 



l-x ' 

Ao stands for the normalization constant, and ySo is the oscillation factor. For the parameters we adopt the 



values from 145], which can describe the data of the nucleon form factors. 
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